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In  recent  years  the  U.  S.  Navy  hag 
' use  increasingly  complex  sonar  systems . The 
^systems  depends  on  a detailed  knowledge  of  % 
acoustic  energy  contained  in  the  underwaterjj 
they  employ.  Once  the  behavior  of  the  accnaj 
then  operational  prediction  in  terms  of  eqal 
uarameters  becomes  possible. 


r^fX  NEL  a great  deal  of  research  ^jfifebeen  devoted  to 

obtalnlngian  objective  tinderstanding  of  the^atur^fvf  the  propaga- 
tion of  acoustic  energy  in  the  convergence  zone  ap A the  surface 
channel  paths*  , As  this  tinderstanding  has  grown  *&re  and  more  effort 
has~fceen‘'SSVCfed  to  the  prediction  of  intensity  fields Since  the 
propagation  of  energy  in  the  convergence  zone  and  surface  channel  i 
paths  is  a complex  phenomena  involving  both  acoustic  and  lenvironmoi- 
tai  parameters,  this  effort  would  not  have  attained  the  success^£p 
has  without  access  to  and  utilization  of  modern  compute^  faciLppes 
and  techniijuea.  CStj,..  t-fptTrT'tfrS'^'ss-es;  LI  ) „ , ~ 


This  presentation  Will  consist  of  three  parts: 
acoustic  intensity  program  for  estimating 


cusslon  of«an  < . . _ _ 

zone  propagation  loss  using  ray  theory;  (2) <a  Alsou 
acoustic  intensity  program  for  estimating  propagatl 
with  surface  channel^  using  normal  mode  theory;  and 
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r -*<  3>fc£vork  at  NEL  -e^  suinarlzlng  information  on  THS 
lables  pertinent  to  these  intensity  computations 
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__|  BAY  THEORY  ACOUSTIC  INTENSITY  PROGRAM.  Th#  input  to  the 

convergence  zone  propagation  loss  program  is  point  temperature  and 
salinity  data  taken  from  the  surface  to  the  bottom  and  the  output  is 
the  propagation  loss  as  a function  of  horizontal  distance  and  depth 
sff  for  a given  source  depth.  I llini  ACCICICTl 
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The  computer  program  complex  consists  of  the  six  sub-programs  listed 
on  slide  1.  The  first  program  computes  the  sound  velocity  from  pres- 
sure, temperature,  and  salinity  after  first  converting  depth  to  pres- 
sure. The  second  program  modifies  the  velocity-depth  data  points  to 
correct  for  the  curvature  of  the  earth.  The  third  program  divides 
the  transformed  profile  Into  layers  and  fits  an  analytic  function  to 
each  layer  such  that  the  velocities  and  their  first  depth  derivatives 
are  continuous  at  the  layer  Interfaces.  The  fourth  program  computes 
ranges,  Intensities,  and  travel  times  using  ray  theory  for  surface- 
reflected,  purely  refracted,  and  bottom  reflected  rays  assuming  a 
flat  bottom.  The  fifth  program  processes  all  the  data  computed  In 
the  fourth  program  for  a given  combination  of  source  and  receiver 
depths.  It  Interpolates,  sums  the  intensities  assuming  random  phase, 
converts  to  decibel  loss,  and  adds  In  attenuation  vith  the  final  out- 
put being  the  propagation  loss  for  100-yard  range  Increments  for  each 
selected  receiver  depth.  It  also  edits  the  data  for  plotting  pur- 
poses. The  sixth  program  plots  enough  points  to  define  the  propaga- 
tion loss  contours.  These  six  programs  consist  of  about  13500  com- 
mands and  constants. 


This  computation  provides  a wealth  of  data  on  propagation 
loss  that  can  be  presented  In  many  ways.  Slide  2 shows  a detailed 
summary  of  the  results  of  one  computation  for  a 6000-foot  source 
depth.  Propagation  loss  as  a function  of  range  and  depth  is  contour- 
ed for  5 db  intervals.  The  degree  of  detail  is  obvious  from  the  com- 
plexity of  the  contouring.  Slide  3 is  a simplified  presentation  of 
the  same  data  prepared  for  a specific  purpose.  Only  one  loss  contour 
Is  shown.  The  shaded  area  shows  where  the  propagation  loss  is  less 
than  100  db  and  the  unshaded  greater  than  100  db.  Additional  pre- 
sentations can  be  prepared  from  the  basic  computations  depending 
upon  the  acoustic  problem  under  consideration.  The  important  point, 
however.  Is  that  the  computation  results  In  a degree  of  objective 
detail  heretofore  unavailable.  This  detail  is  then  available  for 
many  specific  uses. 


The  preceding  discussion  has  been  entirely  theoretical  and 
one  might  reasonably  ask  how  well  the  theoretical  results  agree  with 
experimental  observation.  Slide  4 compares  experimental  and  theor- 
etical results  for  an  acoustic  experiment  conducted  off  San  Diego  In 
January  1955*  The  sound  velocity  profile  was  obtained  from  measure- 
ment of  temperature,  salinity,  and  depth  made  during  the  experiment. 
In  this  experiment  the  source  depth  was  385  feet,  the  receiver  depth 
400  feet,  and  the  frequency  530  cps.  The  irregular  line  connects 
thd  experimental  data  points  and  the  solid  line  shows  the  theoretical 
calculations.  Range  In  nautical  miles  Is  shown  on  the  abscissa  and 
propagation  loss  In  decibels  on  the  ordinate.  The  two  figures  are 
continuous  In  range.  Many  additional  experimental  and  theoretical 


Those  of  you  familiar  with  previous  attempts  to  objectively 
describe  the  acoustic  intensity  field  will  recognize  that  this  com- 
putational procedure  represents  a real  advance  in  understanding  the 
behavior  of  acoustic  energy  being  propagated  in  the  convergence  zone 
path. 


NORMAL  MODE  ACOUSTIC  DETERS ITY  PROGRAM.  The  second  com- 
puter program  complex  that  has  been  developed  at  HEL  is  based  on 
normal  mode,  rather  than  ray  theory,  and  deals  with  the  problem  of 
computing  the  propagation  loss  in  surface  sound  channel  paths.  One 
of  the  important  features  of  normal  mode  theory  is  that  propagation 
loss  can  be  calculated  to  any  range  for  sources  and  receivers  below 
the  layer  as  well  as  in  the  layer.  In  other  words,  there  are  no 
shadow  zones  in  the  computed  Intensity  field  based  on  normal  mode 
theory  as  there  are  in  ray  theory. 


The  model  for  this  computation  is  a bi-linear  duct  consist- 
ing of  one  layer  of  water  having  a linear  sound  volocity  distribution 
overlying  a second  layer  containing  a different  linear  distribution 
of  velocity.  The  inputs  are  depth  and  sound  velocity  gradient  in 
the  first  layer,  gradient  in  the  second  layer,  and  absolute  value 
of  the  surface  sound  velocity.  The  output  is  a plot  of  contoured 
propagation  loss  as  a function  of  range  and  depth  for  a given  source 
depth  and  frequency. 


The  method  is  an  extension  of  the  theory  treated  by  Furry 
and  Marsh.  Early  work  by  these  and  other  investigators  was  limited 
because  of  computational  difficulties  which  required  gross  mathe- 
matical approximations.  These  approximations  limited  the  accuracy 
amd  applicability  of  the  theory  to  a marked  degree.  Utilizing  digi- 
tal  computers  it  is  possible  to  avoid  these  approximations.  This 
problem  Involves  computations  in  the  complex  plane  and  consists  of 
about  6000  commands  and  constants.  Those  of  you  acquainted  with 
normal  mode  theory  fluid  its  applications  are  certainly  well  aware  of 
the  mathematical  complications  associated  with  numerical  applications 
of  the  theory. 


Slide  5 is  a sample  of  the  output  from  this  program.  Con- 
toured plots  come  directly  off  the  on-line  high  speed  printer  and  on 
a CDC  160U  computer  take  about  two  minutes  to  produce.  Range  in 
klloyards  is  shown  on  the  abscissa,  depth  in  feet  on  the  ordinate, 
and  propagation  loss  contoured  for  10  db  intervals.  The  Inputs  for 
this  computation  are  surface  velocity  1*920  feet  per  second,  sound 
channel  depth  100  feet,  sound  channel  gradient  1.8  feet  per  second 
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Again  one  might  ask  how  well  these  theoretical  computations 
agree  with  experimental  observation.  Slide  7 summarizes  results  from 
an  experiment  conducted  off  the  California  coast.  In  this  test  the 
surface  channel  was  300  feet  deep,  the  source  depth  55  feet,  the  re- 
ceiver depth  50  feet,  and  frequency  530  cps.  The  dots  show  the  ex- 
perimental data.  The  solid  line  is  the  theoretical  calculation  based 
on  environmental  factors  measured  at  the  source  ship,  while  the  dash- 
ed line  is  based  on  environmental  factors  measured  at  the  receiving 
ship.  Agreement  between  theory  and  experiment  is  quite  good  as  far 
as  the  general  level  of  propagation  losses  is  concerned.  Note  the 
beat  patterns  that  appear  in  both  the  experimental  and  theoretical 
losses.  These  are  caused  by  the  interaction  between  normal  modes. 


Slide  8 is  another  comparison  from  the  same  data  but  for  a 
l*00-foot  receiver  located  about  100  feet  below  the  surface  channel. 
The  experimental  data  drops  below  noise  between  14  and  15  miles  and 
beyond  17  miles.  At  15  miles  the  experimental  loss  is  some  25  db 
greater  than  for  the  50-foot  receiver  shown  in  the  previous  slide. 
Note  that  normal  mode  theory  predicts  this  large  difference  in  levels 
between  in-layer  and  below-layer  receivers  and  also  predicts  about 
the  proper  decay  with  range.  This  contrasts  markedly  with  ray  theory 
which  predicts  a shadow  zone  beyond  3 miles.  This  ability  to  util- 
ize the  more  exact  normal  mode  theory,  in  place  of  ray  theory,  makes 
it  possible  to  obtain  much  more  realistic  confuted  propagation  loss 
fields. 

1 

ESTIMATION  OF  ENVIRONMENTAL  PARAMETERS.  The  two  computer 
program  complexes  Just  discussed  define  the  environmental  parameters 
for  which  information  is  necessary  in  order  to  obtain  realistic  pro- 
pagation loss  information  for  the  convergence  zone  and  surface  chan- 
nel acoustic  paths. 


In  the  National  Oceanographic  Data  Center  are  archived 
many  millions  of  oceanic  temperature  and  salinity  measurements  made 
over  the  years  for  a variety  of  purposes.  This  data  collection  is 
the  basic  oceanographic  data  collection  for  obtaining  the  required 
input  information.  During  the  past  few  years  the  Oceanometrics 
Group  at  NEL  has  been  attempting  to  develop  techniques  and  methods 
of  summarizing  oceanographic  information  required  by  the  Navy's 
acoustic  program.  To  illustrate  this  work  I would  like  to  review 
the  results  of  two  studies. 


The  first  study  deals  with  the  problem  of  summarizing  sea- 
surface  temperature  observations.  Historically  such  data  have  been 
summarized  by  averaging  over  arbitrary  size  areas  and  time  intervals 
and  the  resulting  averages  plotted  and  subjectively  contoured.  At 
NEL  an  attempt  is  being  made  to  approach  this  problem  from  a more 
objective  point  of  view  utilizing  regression  analysis  concepts.  Al- 
though regression  techniques  were  developed  many  decades  ago  they 
have  rarely  been  used  by  oceanographers  because  of  the  complexity 
inherent  in  developing  realistic  models  and  the  magnitude  of  the 
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arithmetic  task  required  to  evaluate  the  necessary  constants.  With 
the  recent  progress  in  computer  development  the  arithmetic  difficul- 
ties are  being  solved  with  the  result  that  it  is  now  practical  to 
use  complex  models. 

I would  like  to  present  the  results  of  one  such  analysis 
of  sea-surface  temperature  measurements  made  in  the  four  one-degree 
latitude  strips  shown  in  slide  9*  The  measurements  made  in  these 
four  strips  were  treated  as  a single  sample  drawn  from  the  area  30° 
to  U9°N  and  extending  seaward  about  700  miles  and  far  an  l8-month 
time  interval  extending  from  1 April  19^9  to  1 October  1950.  The 
total  number  of  observations  in  each  strip  and  their  distribution  in 
time  and  space  are  also  shown  in  the  slide.  In  the  shaded  areas  one 
to  twenty  observations  were  made  and  in  the  unshaded  areas  no  obser- 
vations were  made. 

Slide  10  shows  the  regression  equation  that  was  fitted  to 
the  observed  data.  Three  main  effects  — latitude,  longitude,  and 
day-of-year  — and  three  interactions  — latitude  by  day,  longitude 
by  day,  and  latitude  by  longitude  — were  considered.  This  resulted 
in  a 22-  variable  equation  with  23  constants  to  be  determined  by 
least  squares. 

i 

Slide  11  summarizes  the  statistical  results  of  the  analysis. 
The  number  of  observations  was  807>  the  per  cent  of  variance  explain- 
ed by  regression  was  85.7  per  cent,  the  multiple  correlation  coeffi- 
cient was  0.93,  and  the  standard  deviation  of  the  observations  about 
regression  was  1 .9°F. 


Slide  12  shows  the  location  in  time  and  space  of  971  tem- 
perature observations  made  in  this  area  during  fiscal  year  1950. 

These  observations  were  not  used  in  obtaining  the  regression  equation 
but  were  used  as  a control  to  see  how  well  the  regression  equation 
could  estimate  independently  observed  sea-surface  temperatures . £ The 
difference  between  the  observed  temperature  and  that  computed-from 
the  regression  equation  was  obtained  and  the  results  are  shown  in 
, the  form  of  a histogram  in  slide  13.  The  standard  deviation  of  the 
differences  was  2.3°F  compared  to  1.9°F  for  the  regression  equation. 
The  close  agreement  certainly  suggests  that  the  two  sets  of  data  are 
drawn  from  the  same  population. 

Slide  lU  shows  a sea-surface  temperature  chart  for  8 Nov- 
ember 1950  contoured  from  temperatures  obtained  from  the  regression 
equation.  It  is  noted  that  the  regression  equation  faithfully  por- 
trays the  main  oceanographic  feature,  namely,  the  tongue  of  cold 
water  that  is  the  result  of  upwelling  occurring  along  the  coast. 

The  second  study  I would  like  to  review  is  concerned  with 


the  shape  of  the  vertical  sound  velocity  profile.  This  study  was 
based  on  an  examination  of  approximately  1000  hydrographic  casts 
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taken  during  the  sunner  and  are  essentially  time  Independent.  Each 
set  of  hydrographic  cast  data  was  converted  Into  sound  velocity  and 
plotted  as  a function  of  depth.  A study  of  these  plotted  data  sug- 
gested three  major  types  which  were  categorized  as  type  A single 
type  B double  thermocllne,  and  type  C double  minimum. 

Slide  15  shows  the  four  single  minimum  sub-types.  Sub- 
type  1 and  sub -type  2 have  a single  minimum  at  about  100  meters  with 
sub-type  1 having  a minimum  velocity  less  than  4800  feet  per  second 
and  sub-type  2 having  a minimum  velocity  greater  than  4800  feet  per 
second.  Sub-type  3 has  a minimum  at  about  500  meters  and  sub-type 
4 a broad,  almost  isovelocity,  minimum  deeper  than  750  meters  and 
sometimes  as  deep  as  1200  meters. 

. 

Slide  16  shows  the  two  double  thermocllne  sub-types.  Both 
sub-types  exhibit  a near  surface  shallow  thermocllne  with  a deeper 
thermocllne  at  about  200  meters  in  the  case  of  sub-type  5 and  at 
about  400  meters  in  the  case  of  sub-type  6. 

Slide  17  shows  the  three  double  minimum  sub-types.  In 
sub-type  7 both  minima  are  shallow  with  the  shallowest  minimum  hav- 
ing the  lowest  velocity,  in  sub-type  8 both  minima  are  deeper  with 
the  deeper  minimum  having  the  lower  velocity,  and  in  sub-type  9 the 
ffl-tirima.  are  still  deeper  with  the  deepest  minimum  having  a much 
lower  velocity  than  the  shallower  minimum. 

After  identifying  these  nine  basic  shapes  each  of  the  1000 
profiles  were  compared  with  the  shapes  and  classified  as  to  sub-type. 
If  a profile  shape  was  different  from  these  nine  basic  shapes  it 
was  classified  as  transitional. 

Next  the  shape  numbers  were  plotted  on  a locator  map  to 
see  whether  or  not  large  volumes  could  be  identified  as  having  the 
shape.  The  results  are  summarized  in  slide  18.  Profiles  taken 
in  the  white  areas  were  of  the  indicated  sub-type.  Gray  areas  con- 
tain transitional  profiles.  There  are  three  major  transitional 
areas  — one  off  the  west  coast  of  North  America,  associated  with 
the  California  Current;  a second  in  the  western  Pacific  associated 
with  the  westward  drift;  and  a third  off  Japan  where  there  is  an 
admixture  of  warm  water  from  the  Kuroshio  Current  and  cold  water 
from  the  Sea  of  Okhotsk.  Thus  this  portion  of  the  North  Pacific 
contains  eight  large  homogeneous  volumes  and  three  transitional 
volumes. 

I 

To  illustrate  what  is  meant  by  the  terms  "transitional 
volume”  and  "homogeneous  volume"  I would  like  to  show  two  sections 
of  profiles.  Slide  19  shows  the  location  of  a section  taken  in  a 
homogeneous  area.  This  600-mlle  section  is  located  in  the  sub-type 
9 volume  with  one  profile  taken  in  the  sub-type  6 volume . Slide  20 
shows  the  individual  profiles.  The  area  to  the  left  of  the  4900 
feet  per  second  value  of  velocity  is  shaded.  The  almost  exact 
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similarity  In  shape  over  distances  of  hundreds  of  miles  is  quite  ob- 
vious. 

Slide  21  shews  the  location  of  a 300-mile  section  located 
In  the  transitional  volume  off  the  coast  of  Japan . The  upper  portion 
of  slide  22  shews  the  profiles  over  northern  half  of  the  section. 
These  shapes  exhibit  relatively  small  shape  changes.  The  lower  por- 
tion shows  the  profiles  for  the  southern  half  of  the  section  and 
Illustrates  marked  changes  in  profile  shape  occurring  over  distances 
of  20  to  bO  miles  — approximate  convergence  zone  distances. 

This  analysis  suggests  three  Implications  --  to  operations 
It  implies  that  In  the  homogeneous  areas  acoustic  performance  should 
be  stable  and  predictable  while  in  the  transition  areas  It  should  be 
variable  and  less  predictable;  to  research  it  suggests  that  acoustic 
experiments  should  be  conducted  in  the  different  profile  volumes  and 
the  transitional  volumes  before  generalizations  are  attempted;  and 
to  survey  programs  such  analyses  gives  guidance  as  to  areas  that 
should  be  surveyed. 

PREDICTION  OF  ACOUSTIC  PERFORMANCE.  In  the  discussion  so 
far  Information  explaining  acoustic  results  by  means  of  theoretical 
calculations  based  on  oceanographic  data  taken  during  acoustic  tests 
in  familiar  waters  has  been  presented.  Our  current  ability  to  pre- 
dict the  acoustic  field  In  areas  significantly  different  from  NEL's 
customary  test  areas  is  suggested  by  slide  23.  Convergence  zone 
echo-ranging  experiments  were  conducted  at  two  locations  --  area  1, 
located  in  a transition  region  300  miles  west  of  Vancouver  Island 
and,  area  2,  located  in  the  homogeneous  region  100  miles  south  of 
Kodiak.  In  laboratory  preparation  for  these  tests,  appropriate 
oceanographic  stations  were  examined  and  sound  velocity  profiles 
forecast  for  the  two  areas.  Predicted  propagation  losses  were  cal- 
culated for  each  of  these  areas.  Three  horizontal  bars  are  shown 
for  each  area.  The  upper  bar  shows  the  range  interval  at  the  first 
convergence  zone  over  which  echoes  were  experimentally  observed,  the 
next  bar  the  expected  Interval  calculated  by  using  a sound  velocity 
profile  obtained  from  measurements  made  concurrent  with  the  accxistle 
experiments,  and  the  lower  bar  the  predicted  interval  obtained  from 
a predicted  sound  velocity  profile.  Range  prediction  obtained  from 
the  predicted  sound  velocity  profiles  was  not  as  good  in  the  trans- 
itional volume,  area  1,  as  In  the  homogeneous  volume,  area  2,  pre- 
sumably because  of  large  horizontal  variation  in  profile  shape. 
However,  in  both  Instances  the  agreement  between  the  predicted  detec- 
tion range  and  the  observed  detection  range  Is  quite  remarkable  at- 
testing to  the  validity  of  the  prediction  approach. 

I would  like  to  conclude  with  a comment  on  the  possibility 
of  developing  an  operational  acoustic  performance  prediction  tech- 
nique based  upon  the  research  I have  Just  discussed.  From  what  has 
been  accomplished  so  far  it  seems  well  within  the  realm  of  possibil- 
ity to  develop  a system  for  predicting  the  sound  velocity  distribu- 
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tion  and  its  Influence  on  various  acoustic  parameters,  from  periods 
of  a few  days  to  a few  tens  of  days;  as  well  as  to  give  the  oceano- 
metrlc  expectations  for  any  time  period  in  the  future.  It  is  not 
very  difficult  to  visualize  a shipboard  equipment  where  the  analyzed 
historical  oceanographic  data  and  sonar  equipment  parameters  have 
been  stored;  observations  taken  during  the  preceding  few  days  or 
weeks  introduced;  and  by  entering  the  latitude,  longitude,  and 
time-of-year  obtaining  a predicted  acoustic  field. 
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RAY  THEORY  ACOUSTIC  INTENSITY 
COMPUTATIONAL.  PROGRAM 


Sub- Program 

1.  Sound  Velocity 

2.  Curved  Earth  Simulation 

3.  Curve  Fitting 

4.  Acoustic  Intensity 

5.  Propagation  Loss  Computation 

6.  Propagation  Loss  Plotting 
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RECEIVER  DEPTH:  400  FT. 
FREQUENCY:  SSOCPS. 


RANGE  (NAUTICAL  MILES) 


Figure  8.  Experimentally  Observed  Propagation  Loss  Compared  with 
Propagation  Loss  Computed  from  Normal  Mode  theory 
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MULTIPLE  REGRESSION  EQUATION 
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STATISTICAL  RESULTS  OF  THE 
REGRESSION  ANALYSIS 


NUMBER  OF  OBSERVATIONS  ®07 

PERCENT  VARIANCE 

EXPLAINED  BY  REGRESSION  85.7* 

MULTIPLE  CORRELATION  COEFFICIENT  0.83 

STANDARD  DEVIATION  OF  THE 

OBSERVATIONS  ABOUT  REGRESSION  1.9  F 


Figure  11. 
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Used  to  Verify  Regression 
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Figure  20.  Sound  Velocity  Profiles  Taken  in  a Homogeneous  Volume 
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Figure  21.  Location  of  Sound  Velocity  Profile  Section 
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Figure  22.  Sound  Velocity  Profiles  Taken  in  a Transition- Voluae 
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Figure  23.  Experimental  Verification  of  an  Acouatic  Range  Prediction 


